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Abstract 

We construct static black hole solutions that have no rotational symmetry. These arise 
in theories, including the standard electroweak model, that include charged vector mesons 
with mass m 7^ 0. In such theories, a magnetically charged Reissner- Nordstrom black hole 
with horizon radius less than a critical value of the order of is classically unstable 
against the development of a nonzero vector meson field just outside the horizon, indicating 
the existence of static black hole solutions with vector meson hair. For the case of unit 
magnetic charge, spherically symmetric solutions of this type have previously been studied. 
For other values of the magnetic charge, general arguments show that any new solution with 
hair cannot be spherically symmetric. In this paper we develop and apply a perturbative 
scheme (which may have applicability in other contexts) for constructing such solutions in 
the case where the Reissner-Nordstrom solution is just barely unstable. For a few low values 
of the magnetic charge the black holes retain a rotational symmetry about a single axis, 
but this axial symmetry disappears for higher charges. While the vector meson fields vanish 
exponentially fast at distances greater than 0(m^^), the magnetic field and the metric have 
higher multipole components that decrease only as powers of the distance from the black 
hole. 
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1. Introduction 



One of the many remarkable features of black holes is the symmetry and simplicity of 
the time-independent black hole solutions. The static vacuum black holes are all spherically 
symmetric and determined by a single parameter. Adding electromagnetism gives the pos- 
sibility of endowing the black hole with electric or magnetic charge, but the static solutions 
remain spherically symmetric, with purely Coulomb electromagnetic fields. Even if one con- 
siders solutions that are stationary, but not static, a rotational symmetry about one axis 
remains. This situation stands in constrast with that of electromagnetism in flat spacetime, 
which possesses static (even if singular) solutions corresponding to point multipole moments 
of arbitrarily high order. An explanation for the absence of static gravitational solutions with 
higher multipoles comes from the no-hair theorems [1] that sharply constrain the possible 
structure of black holes both in the electrovac case and for gravity coupled to a number of 
types of matter. 

However, it has become clear in recent years that if the theory governing the matter fields 
has sufficient structure, it is in fact possible to have black holes with nontrivial static fields 
outside the horizon; i.e., black holes with hair. In particular, theories with electrically charged 
massive vector mesons can have two types of magnetically charged black hole solutions [2, 3] . 
One is the trivial generalization of the Rcissner-Nordstrom solution to the coupled Einstein- 
Maxwell equations. The other, which exists only if the horizon radius is sufficiently small, 
has nonzero massive vector fields just outside the horizon. For the case of the SU{2) gauge 
theory with a triplet Higgs field, which has a nonsingular magnetic monopole solution in 
flat spacetime, one finds a new solution with unit magnetic charge that may be viewed as 
a Schwarzschild-hke black hole embedded in the center of an 't Hooft-Polyakov monopole 
[4] . Although this solution was first found directly, a signal of its existence is the fact that 
the Reissner-Nordstrom black hole develops a classical instability when its horizon radius 
becomes smaller than the radius of a magnetic monopole core [5] . 

Similar arguments [6] based on instabilities of Reissner-Nordstrom solutions suggest the 
existence of new black hole solutions with higher magnetic charges. However, in the presence 
of a magnetic monopole a spherically symmetric charged spin-one field is possible[7,8] only 
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if the product of the magnetic charge of the monopole and the electric charge of the field 
is unity Hence, these new black holes can be at most axially symmetric. Whether or not 
they actually possess such symmetry is a question not of general principle, but of detailed 
dynamics. In this paper we will show that, at least for certain ranges of parameters, they 
do not. 

A theory with sufficient structure to yield these new black holes has matter fields de- 
scribed by the fiat-spacetime Lagrangian 

where 

F^,^. = d^,A^ - d^^A^ (1.2) 



W^^ = D^W^ - D^W^ (1.3) 

D^,W,^{d^,-ieA^)W,. (1.4) 

The fourth term in the Lagrangian is an anomalous magnetic moment term, with the constant 
g arbitrary. In order that the energy be bounded from below, we must require that A > / 4 
[6]. 

If we add to the theory a neutral scalar field with appropriate self-interactions and 
give the vector field a ^-dependent mass m — ecp, then for = 2 and A = 1, then Eq. (1.1) is 
simply the unitary gauge form of the Lagrangian for an SU (2) gauge theory spontaneously 
broken to U(l) by a triple Higgs field. Similarly, for = 2, A = 1/ sin^ 6w, and m — e0/2 we 

1 This follows from the absence of vector spherical harmonics with total angular momentum zero. The 
absence of such harmonics can be understood by considering the motion of a particle with electric 
charge e about a monopole with magnetic charge Qm- The total angular momentum is the sum of 
the spin angular momentum, the orbital angular momentum, and a contribution of magnitude cQm 
directed along the line from the monopole to the charge; if cQm 1) the sum of these three terms can 
never vanish. 
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obtain the unitary gauge form of the standard electroweak Lagrangian, but with all terms 
involving the Z or fermions omitted. It is a straightforward matter to extend the analysis 
of this paper to such models. 

We are seeking static black hole solutions to the theory obtained by coupling the La- 
grangian of Eq. (1.1) to general relativity. In the absence of rotational symmetry, the static 
field equations are a set of coupled partial differential equations in three variables. An exact 
analytic solution of these is beyond our abilities. Instead, we use a perturbative approach. 
We begin by considering a Reissner-Nordstrom black hole with radial magnetic field 

Fe<^ = ^sin^ (1.5) 

corresponding to a magnetic charge Qm = q/g (where q is restricted by the Dirac quantiza- 
tion condition to integer or half-integer values) and vanishing W field. The metric is 

ds'^ = -B{r)dt^ + B-\r)dr'^ + r^{de^ + sin^ Odcj)'^) (1.6) 

with 

AM = l-gj£g + j^ =('■-■«)(''-''-). (1.7) 

We choose the mass M to be such that the outer horizon radius ru is less than rcr, the 
critical value for instability. It is in this mass range that we expect there to be a second black 
hole solution with nontrivial W field and with electromagnetic field strengths and metric that 
differ from the Reissner-Nordstrom form. It is often the case that the exponentially growing 
eigenmodes about an unstable static solution give a good indication of the nature of a nearby 
stable solution, particularly in the case where the original solution is just barely unstable. 
Guided by this intuition, we linearize the static field equations about the Reissner-Nordstrom 
solution. This leads to an eigenvalue problem that is closely related to, although not exactly 
the same as, that encountered in the stability analysis [9]. For ru close to rcr, there is 
a single negative eigenvalue, whose magnitude tends to zero as r^/ — > rcr. With rcr — tu 
sufficiently small, this eigenvalue becomes a small parameter that can serve as the basis for 
a perturbative expansion. 
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In Sec. 2, we illustrate our method with a simple toy model consisting of a scalar field 
coupled to a fixed, but spatially inhomogneous, external source. In Sec. 3, we set up the 
formalism for treating the case in which are actually interested, that of a black hole in 
the theory described by the Lagrangian of Eq. (1.1). We assume that Gm^/e^ is small; this 
allows one to solve (to leading order) for the charged vector field and the perturbations of the 
electromagnetic field before dealing with the metric perturbations. For technical reasons, it 
turns out that the details of the subsequent analysis are considerably simpler if g is positive 
and q > 1. We exploit these simplifications in Sec. 4, where we determine the leading 
perturbations of the electromagnetic field in terms of those of W^. In Sec. 5, we examine 
the lowest order contributions to the charged vector field and show that there is a parameter 
range for which the solution is not even axially symmetric. In Sec. 6 we obtain the leading 
corrections to the Reissner-Nordstrom metric. Some conchiding remarks are included in 
Sec. 7. An appendix describes some results needed to construct Green's functions that we 
use. 

2. A Toy Model 

Consider a real scalar field whose dynamics is governed by the Lagrangian 

^--l {d,^f - ^^(x)0^ - \ct>^ (2.1) 

where F(x) arises from the coupling to a spatially inhomogeneous but static external source. 
Static solutions obey 

= [-V^ + F(x)]0 + X<f 

(2.2) 

and, to have finite energy, must satisfy the boundary condition that vanish at spatial 
infinity. 

The trivial configuration 0(x) = is a solution for any choice of F(x). It is easy 
to see that this is the only static solution if Al is a positive operator. However, if M. 
has negative eigenvalues, this trivial solution is unstable, implying the existence of a new. 
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spatially inhomogeneous, solution whose form we seek. To this end, let us assume that A4 
has only a single negative eigenvalue, with an eigenf unction tjj obeying 



and normalized so that 



(2.3) 



(2.4) 



We now write as the sum of a term proportional to ip and one orthogonal to it; i.e.. 



with 



(f){x) — kip{x) + 0(x) 



J d^xip{x)4>{:x.) = 0. 



The static field equation (2.2) then implies that 
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M(j) + X[kip + (j)) +Ti! = 



(2.5) 



(2.6) 



(2.7) 



where F is a Legendre multiplier that enforces the orthogonality condition (2.6). It can be 
calculated by multiplying both sides of this equation by ip and then integrating over all space 
to obtain 



r = -A / d-^x 



Variation of the action with respect to k gives the additional equation 

91 . 
dk-' 



(2.8) 



(2.9) 



where 



7 = j d^x ^k'^ipMip + ^ (^kij + (f) 



(2.10) 



Thus far we have made no approximations; together, Eqs. (2.7), (2.8), and (2.9) are 
completely equivalent to Eq. (2.2). We now recall that if 6 = the scalar field 0, and 
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therefore k and 0, must vanish. Hence, for small b it should be possible to expand these 
quantities as power series in b. Furthermore, since it is the existence of the negative eigenvalue 
which makes a nontrivial solution possible, we may view kij) as providing the source for 
(through Eq. (2.7)). We therefore expect that cp is of higher order in b than kip. Assuming 
this to be the case, Eqs. (2.9) and (2.10) give 



I 



+ 0{b^) 



(2.11) 



while Eq. (2.8) implies 



r^-Xk^ J d^x'ilj^{^)+0{b'^) 



d^xijj'^ix) 



.-1/2 



+ Oib') . 



These may be substituted into Eq. (2.7) to give 



(2.12) 



M4>{x) = -Xk" 

~ Vx 



V^3(x)-V^(x) J d'yip\y) 

■3/2 



d^xiP'^ix) 



+ 0{b^) 
V^3(x)-V^(x) J d^yip\y) 



(2.13) 



+ 0{b^). 



This shows that is of order b^, and justifies the assumption above that it is of higher order 
than kip. Thus, to leading order the static solution is approximated by the negative eigen- 
value fluctuation about the vacuum solution, multiplied by a scale factor whose magnitude 
is determined by the nonlinear term in the Lagrangian. 

By substitution of the lower order results back into the original equations, k and 0, and 
hence itself, can be calculated to to arbitrarily high order in b. 
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3. Charged Vector Meson Model 



We now apply this method to the theory in which we are actually interested, namely 
that described by the Lagrangian of Eq. (1.1). The first step is to identify the unstable (i.e., 
exponentially growing in time) modes about the unperturbed Reissner-Nordstrom solution. 
This was done in Ref. 9, whose results we now briefly summarize. When the field equations 
are linearized, the perturbations in the gauge field and the metric decouple from those in 
the massive vector field. The linear perturbation problem for the former two modes is the 
same as in the pure Einstein-Maxwell theory, where it was shown some time ago [10] that 
the Reissner-Nordstrom solution is stable. Hence, the stability analysis reduces to a study 
of the linearized W field equations. It is convenient to define M-n^ by 

M^'^W,, = -^Da (x/IW^"^) + m^W" - '-^F^'^'Wa (3.1) 

where here, and for the remainder of the paper, we adopt the convention that g^iy, A^, 
and Fij_i, denote the corresponding unperturbed quantities while is the gauge covariant 
derivative taken with respect to the unperturbed potential; indices are raised and lowered 
with the unperturbed metric. The unstable modes are solutions of 

M^'W^ = (3.2) 

whose time-dependence is of the form 

Py^(x,t)=/4x)e'^* (3.3) 

with real lo. The spherical symmetry of the unperturbed solution allows one to choose the 
solutions of Eq. (3.2) to be eigenfunctions of both and J^, where J is the total angular 
momentum operator. Because of the extra angular momentum of an electric charge in the 
field of a magnetic monopole, the corresponding eigenvalues are not the usual ones. Instead, 
J runs in integer steps upward from the minimum value Jmin = ? — 1, unless q — 1/2, in 
which case Jmin = 1/2. For each value of J, unstable modes exist if the horizon radius rjj is 
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less than a critical value rcr(J) that is of order m~ , provided that g lies in an appropriate 
range {g > foi J = q — 1, g > 2 foi J = q, and either g < or g > 2 for J > q). For 
a given value of g, rcr(J) is greatest for the smallest J that can have unstable modes with 
that g. Thus, if J denotes the value that maximizes rcr, we have J — q — lfoTg>0 and 
J = ? + 1 for ^ < for 5 > 1. If 5 = 1/2, J = 1/2 if > 2 and 3/2 if ^ < 0; if < ^ < 2, 
there is no instability. 

The modes that will form the basis for our new solutions are the static eigenfunctions of 
M. with negative eigenvalue; i.e., the time-independent solutions of 

M/V'. = -/3V^|J^ (3.4) 

with real /3. (A factor of has been extracted to make /3 dimensionless.) 

This eigenvalue equation must be supplemented by boundary conditions. At spatial 
infinity we merely require that ip^ not diverge. For negative eigenvalues (indeed for all 
eigenvalues less than m^) this implies that in fact vanishes as r ^ oo. A second boundary 
condition is obtained at the horizon, where we require that ip^ be regular, in the sense that its 
components measured relative to a coordinate system that is nonsingular at the horizon (e.g., 
Kruskal-like coordinates) be regular. Because of the manner in which the singular metric 
factors enter Eq. (3.4), this constrains the behavior of ip^ near the horizon — as we will see 
more explicitly in the next section — and causes the spectrum of negative eigenvalues to be 

2 

discrete. For this portion of the spectrum, we can require that the eigenfunctions satisfy the 

3 

normalization condition 

Jd^xv§r^,r-i (3.5) 

where, both in this equation and hereafter, the spatial integration is understood to be re- 
stricted to the region outside the Reissner-Nordstrom outer horizon. 

2 It might seem strange that the nature of the spectrum should be determined by the singularities of a 
metric at a horizon that is only a coordinate singularity. This happens because the condition we are 
imposing on the eigenfunctions, that they be static, is defined in terms of a coordinate t that is singular 
at the horizon. 

3 Since, as is easily shown, static solutions of Eq. (3.4) must have V't = 0, tp^ip^ is positive. 
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Because of both the nontrival metric component gtt{r) and the possibihty of a nonvan- 
ishing Wt in the time-dependent case, the eigenfunctions are not in general the same as 
the that appear in Eq. (3.3): the spectra of the uj and (3 are not even the same. How- 
ever, a zero eigenvalue for the static operator does correspond to a zero frequency of the 
small oscillation problem and, furthermore, the static problem has negative eigenvalues if 
and only if the Reissner- Nordstrom solution is unstable. Hence, the conditions for instability 
enumerated above are also the conditions we need to be able to construct our new solutions. 

In fact, the static eigenmodes for real (3 can be obtained from the uj = solutions of 
Eq. (3.2) with different values for the parameters. This can be seen by bringing the right 
hand side of Eq. (3.4) over to the left; the resulting equation is precisely that satisfied by 
/^(x) for a; = 0, but with replaced by m^(l + /3^). It follows that the value of rn that 
leads to a given /3 for l^-mass m is equal to the critical value rcr for a VF-mass my^l + f3^. 
If m <^ Mpi (the case with which we will be primarily concerned), rcr is much greater than 
the horizon size for an extremal Reissner-Nordstrom black hole and its dependence on the 
inner horizon r_ can be neglected. Dimensional arguments then show that rcr is inversely 
proportional to m. It then follows that 



P 



— r 



H 



(3.6) 



Thus, by taking ru — rcr (J) <S fcx{j)-i we ensure that /3 <^ 1, thus providing the small 
parameter needed for our perturbative calculation. 

In general J is nonzero, so that instead of a single unstable mode, as in the model of 
Sec. 2, there is a degenerate multiplet of unstable modes -0^^ that are distinguished by the 
eigenvalue of J^. Proceeding as in that section, we write the W field as a linear combination 
of the unstable modes plus a remainder orthogonal to these modes, 

J 

W^,^Vf, + W^,^m-^l'' + (3.7) 

M=-J 
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where 

J d'xV§r^W^-0. (3.8) 

It is useful to define a quantity a by 

J 

E \kM?^a' (3.9) 

M=-J 

SO that 

= 0(a) . (3.10) 

Since, will be displayed explicitly below, the source for the perturbations of the electromag- 
netic field is quadratic in and contains an explicit factor of e, 

SAi^ = 0(ea2) . (3.11) 

(Note that the quantities ea and Gm^/e^ are truly dimensionless, whereas e is dimensionless 
only if one sets h = 1, which would not be natural in this essentially classical context.) 
The source for the metric perturbations is the perturbation of the energy-momentum tensor. 
The leading contribution to this, of order a^, is from terms quadratic in and from terms 
linear in SAfj,. However, these enter the field equation multiplied by a factor of G, and so 
= /i^jy must be suppressed by an additional factor of roughly Gm^ — (m/Mpi)^. Hence, 

= V = 0{Gm'^a^) . (3.12) 
Finally, the magnitude of W can be determined from the field equation 

ie ^ ^ (3.13) 

+ — {V^bA- - V^bA% - YmVm + • • • ■ 

Here the dots represent terms which are either O(e^a^) or 0{Gvr?c?^ or smaller, while the 
V M are Lagrange multipliers introduced to enforce the orthogonality of W^j, and the -0^. One 
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can solve for the Tm by multiplying both sides of this equation by and then integrating 
over all space outside the horizon. Inserting the result back into Eq. (3.13), we see that 



= 0{e'^a^) . (3.14) 

We will assume that Gm? /e^ <^ 1. The leading behavior of V^, VF^, and 5A^ can then 
be obtaining by solving the field equations in the background of the unperturbed Reissner- 
Nordstrom metric. Having done this, the leading perturbations of the metric can then be 
obtained. In fact, for calculating the lowest order metric perturbations, only dA^ and 
are needed. For the remainder of this section, and the next two, we will concentrate on the 
determination of these two quantities. We will then return to the calculation of the metric 
perturbations in Sec. 6. 

Linearization of the electromagnetic field equation about the unperturbed solution yields 

^d, {^5F^^) = -^d, (V^p^-) +f (3.15) 



where 



P,u^--^{V:V,-V:v,) + .-- (3.16) 



and 



f = ie [V* {D^V - D^V) - Vf, {D^V*" - D^V*!")] + • • • (3.17) 
with the dots signifying higher order terms. Similarly, the Bianchi identity gives 

e^"'"'f^dJFai3^0. (3.18) 



In addition to these, we need the equations obtained by varying the action with respect 
to the kM- In deriving these, we need only take into account terms in the action of up to 
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order e a and so can approximate the action by 



approx — 



V 



(3.19) 



(Terms linear in both and Vy^ which would be of order e^a^, are absent because of the 
orthogonality of and the ^ There are two ways of proceeding from here. One can 
substitute the expansion of Eq. (3.7) for and then vary the above expression with respect 
to /cm, thus obtaining an equation involving both the and bAy. Alternatively, one can 
first use Eqs. (3.15) and (3.18) to solve for bA^^ and bF^^^y in terms of and then substitute 
the resulting expressions back into Eq. (3.19) to obtain an action which is a function of only 
the we will follow this second approach. A number of simplifications arc possible. First, 
the eigenvalue equation (3.4) and the normalization condition (3.5) can be used to integrate 
the term quadratic in T/^. Next, by multiplying both sides of Eq. (3.15) by ^bAy and then 
integrating by parts, one obtains the identity 



which can be used to eliminate the term quadratic in 8F^y. A similar procedure applied to 
the source-free equation obeyed by the unperturbed field strength shows that 



(In both cases, one can verify that the surface terms from the integration by parts vanish 
as long as the total magnetic charge is held fixed.) The term quadratic in Fy^y is obviously 
independent of the ku and can be ignored. Finally, since all quantities are independent of 
time, it is sufficient to integrate over the spatial variables. We are thus led to the equation 




(3.20) 




(3.21) 



dl 



(3.22) 



= 



dku 



12 



where 



Ae2 , 



(3.23) 



and 5 Ay is understood to be given in terms of the kM- Since the integrand on the right hand 
side is of order e^a^, we see that a is proportional to /3/e, indicating that our perturbative 
expansion is justified for ru sufficiently close to rcr- 



4. The Case q > l,g > ^ 

We now speciahze to the case g > 1, 51 > for which, as was noted above, J — q — 1. 
This allows us to take advantage of the special properties [8] of the J — q — 1 vector spherical 
harmonics, which lead to a number of technical simplifications in the analysis. For J = q — 1, 

and only for that value, there is but a single monopole vector spherical harmonic [8,11] for 
each value of = M. Hence, if we denote this harmonic by Cj^{9,(p), the unstable modes 
in this case can be written in the form 

V'f = /(r) (4.1) 

where /(r) does not depend on M. 

The J — q — 1 harmonics have a number of special properties. Their radial and time 
components vanish, 

C^^ = (7/^ = 0, (4.2) 
and their two angular components are related by 

Cf ^isinOC^ . (4.3) 

In addition, their covariant curl, evaluated in the background Dirac vector potential, van- 
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ishes: 



M 



as does their covariant divergence 



1 



0. 



(4.4) 



(4.5) 



A convenient choice of normahzation condition is 



/ 



dcf) d9 sin 6 



Cf(^,0)l C^^(e,0) 



(4.6) 



To obtain an exphcit expression for the , we must choose a gauge. If the electromag- 
netic vector potential has a single nonvanishing component 



A<. = ^(l-cos^), 



then 



where 



sin 9 



1 + cos 61 



2'?V27r 



(2g - i; 



q+M-l 



1/2 



Xq + M-l)\{q-M 
With the aid of these properties, the eigenvalue equation (3.4) reduces to 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



where /(r) can be chosen to be real. Equations (3.5) and (4.6) fix the normalization of / to 
be 



(ir|/(r)|2=l. 



(4.11) 



Given rcr, and hence /3, Eq. (4.10) can be integrated numerically to obtain /(r). This 
function is monotonic, has no zeros, and vanishes exponentially with r as r — > oo. Near the 
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horizon it behaves as 

/(r) = A[l - b{r - rn)] + 0[{r - rnf] 



(4.12) 



where 

[B'{rH)r'>0. (4.13) 



19 2/1 , n2 

-2 m (1 + p 



Proceeding with the construction of the solution, we write 

V^^m-^/''f{r)^^{eA) (4.14) 

where 

$^(^,0)= kMCff{e,<t>). (4.15) 

M=-{q-l) 

Eqs. (4.2) and (4.3) imply that $.r = $t = and fix the ratio of $61 and $0. Note that <l>^ 
has exactly 2(g — 1) zeros as 6 and range over the unit sphere. To show this, we use the 
explicit expression (4.8) for the vector harmonics and write 

q-l 

^Q{e,(t>) = e''^{l + coseY-^ aqMkMz'^+^-^ (4.16) 

M=-{q-l) 

where 

smS (4.17) 
1 + cos 61 ^ ^ 

The entire complex 2;-plane maps onto the unit sphere, with |z| = oo corresponding to the 
south pole, 6 = 71. Let M be the largest value of M for which kM is nonzero. The sum in 
Eq. (4.16) is then a polynomial of order M + q — 1 in z, and thus has M + q — 1 zeros at 
finite z. In addition, the pref actor multiplying the sum combines with the M — M term to 
give a zero of order g — 1 — Mat^ = 7r. Adding these together, we obtain the promised 
result. 
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The properties of the Cjf also lead to the useful identity 

- K% = ir%.K^'' (4.18) 
where e^jy is an antisymmetric tensor whose whose only nonzero components are 

= sin^. (4.19) 



We will encounter the quantity in the source terms for the perturbations of both 

the electromagnetic field and of the metric. Since we will solve these equations by separation 
of variables, it is useful to define the expansion 



r2 



$;(^,0)$^(^,0) = a^Y.''j^^irn{0,<P) (4.20) 



jm 



where 

ajm = d<pd0 sine Y*^{e,<P)%{e,<f>)^^'{0,<f>) (4.21) 

and a is as defined in Eq. (3.9). Ordinary, rather than monopole, spherical harmonics enter 
here because we are dealing with a neutral quantity. Hence, j runs over integer values 
although, since is a linear combination of monopole harmonics with angular momentum 
J — q—1, the ajm vanish for all j > 2(^—1). Note that, as a consequence of the normalization 
condition (4.6), 

(TOO = ^ . (4.22) 
v47r 



The properties of the J — q — 1 harmonics also simplify the electromagnetic field equa- 
tions. All components of vanish, while 

P,. = ^^M'^^V'^;^" ■ (4.23) 



16 



The various components of Eq. (3.15) can be written as 



d^{^dF''')=0 (4.24) 

9^(V?<^Fn=0 (4.25) 

-^d, iV^SF^^-) = ($*$") (4.26) 

where we have adopted the convention that Roman indices from the beginning of the alphabet 
take only the values 9 or 0. (In obtaining the last of these, we have used the fact that ^/g e°'^ 
is a function only of r.) 

Because wc arc seeking time-independent solutions, the equations for the electric and 
magnetic fields decouple. For the former the source term vanishes, and so the equations are 
the same as those encountered in studying perturbations of the pure Reissner-Nordstrom 
solution, where the only allowed static perturbation of the electric field is a radial field 
corresponding to a variation of the black hole's electric charge. Since we are assuming 
vanishing electric charge, this perturbation must be excluded, and so = 0. 

The equations for the magnetic field can be solved by separation of variables. We first 
expand SF^p in terms of vector spherical harmonics: 

5Fe^ = Y,Fl'^{r) smeYj^{e,<t>) 

r 1 (4-27) 

6Fra = J2 [Frir)ea%Yjm{e, 0) + F^ir) daYjm{e, 0) 

jm 

where it is understood that = F^^ — 0. Substituting this expansion into Eq. (4.25) and 
using the identity 

^ {V^r^d.Yjm) = -^-^^-i^Yjm (4.28) 

we find that 

fI"^ = . (4.29) 
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Next, the /i = t component of the Bianchi identity, Eq. (3.18), leads to 



^ = iil±llFr . (4.30) 
dr ^ 

For j > 1 this can be used to ehminate F^^, while for j = it implies that F^^ is a constant. 
Since a constant corresponds to a change in the magnetic charge, we set = 0. 



Finally, Eq. (4.26), together with Eqs. (4.20) and (4.30), yields 

X B—J:- Z^^i = -T^-^U + 1)/ ^jm ■ (4.31) 



It will be convenient to write 



where Tj obeys 



f-rW = ^<T^„^:/(r) , J>0, (4.32) 



^U'!£i].i(l±}lr, = -i(i^f. (4.33) 
dr \ dr mr^ 



By multiplying this equation by Tj and then integrating over r from to oo, one can 
show that for / = and j > the only regular solution is the trivial one J^j{r) = 0. Hence, 
in the presence of the source one can solve, at least formally, for J-'j by inverting the operator 
on the left hand side of Eq. (4.33). To construct the appropriate Green's function we need 
the two solutions gj'{r) and gj~{r) of the homogeneous equation that are regular at r — rjj 
and r — oo, respectively. Using the fact that B{r) tends to unity at large r, one immediately 
finds that these two solutions behave asymptotically as r^'^^ and r~^. (Explicit forms for 
these solutions are given in the appendix.) If they are normalized so that r~^^'^^^g~{r) and 



r^g~j{r) both tend to unity as r — > oo, then 



00 

TAr)^-'-^jdr'Gf\ry) 



(4.34) 



rH 
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where 



Gf\r, r') = - [e{r - r')g+{r)g-{r') + e{r' - r)g-{r)g+{r')\ . (4.35) 

Note that neither g^{r) nor gj{r) can have any zeros for r > rjj (i.e., in the region where 
B{r) > 0). This fact, together with Eqs. (4.34) and (4.35), imphes that J-j{r) is positive 
everywhere outside the horizon. 

Because /(r) falls exponentially for r ^ m~^, the contribution from the first term in 
the Green's function dominates at large distance and so 

JT^- ~ ^ , r ^ oo (4.36) 



rJ 



where 



oo 



'J 



2j + 

rH 



tl / dr ^ ^ ' . (4.37) 
1 7 mr2 ^ ' 



As a check that this is indeed the proper behavior, note that this imphes that the 2-^-pole 
components of 5Fq^ fall as 1/r-^, while the unperturbed monopole component of is 
independent of r. Since a is of order /3/e, the magnetic field perturbations that we have 
found correspond to magnetic 2-^ -poles with components equal to the ajm times quantities 
of order ~ 0^ jem? . 



5. Determination of the ku and the symmetry of the solution 

We can now use our results for and bF^^^y to determine the ky[. The overall scale of 
these, measured by the quantity a that was defined in Eq. (3.9), determines the magnitude 
of the departure of our solution from the Reissner-Nordstrom black hole. The relative sizes 
of the various ku determine the angular dependence — i.e., the shape — of the solution; for 
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studying these it is convenient to define 

riM^— (5.1) 
a 

that satisfy 

^|nMp = l. (5.2) 

M 

Substitution of the of the results of the previous section into Eq. (3.23) gives the quantity 
I ^-(3 ma + p 2^ \ajm\ ^ 2^ Qj 2^ Wjm\ (5.3) 

jm j m 

whose minimum determines the kM- In this expression p and qj denote the positive integrals 
and 

= / ^^{^^^W (5.5) 

which are all of order unity. 

Eqs. (4.15) and (4.20) show that the ajm are homogeneous polynomials of degree 2 in 
the Um- This suggests that we rewrite Eq. (5.3) as 

/ = -/3W + ^^4^/i(nM) (5.6) 

o 

where 

2(g-l) 2(g-l) 

Hum) = ^P (^•'^) 

and the rotational scalars 

j 
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are homogeneous polynomials of degree 4 in the um- Minimization of / requires 



a = ^ [hiriM)]-'^^ . (5.9) 

eg 



The Um are determined, up to an ambiguity corresponding to the rotational and global 
gauge symmetries of the theory, by minimizing Ii. We begin by considering individually 
several low values of q. 

Case i: q = 1, J = 

The solution with unit magnetic charge is spherically symmetric (indeed, it is the only 
case for which spherical symmetry is possible). There is only a single um, of unit magnitude, 
whose phase has no physical significance. 

Case ii: q = 3/2, J = 1/2 

There are two um that form a complex SU{2) x f/(l) doublet, where the former factor 
refers to spatial rotations and the latter to global phase rotations of the charged fields. It is 
always possible to find a symmetry transformation that brings such a doublet into the stan- 
dard form 774/2 = Ij "'-1/2 = 0- The solution is axially symmetric in the sense that it is left 
invariant by a combination of a rotation about the z-axis and a global gauge transformation. 
In particular, all gauge-invariant quantities are manifestly axially symmetric. One finds that 

r2$*$M = _(i_cos^) (5.10) 



so that the nonzero Ujm are 



(700 = ^=, (710 = 7==- (5.11) 



The solution has a net magnetic dipole moment that can be attributed to the asymmetric 
distribution of the magnetic dipole density of the charged vector field. 
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Case iii: q = 2, J = 1 

This case is somewhat less trivial. The three complex um are equivalent to a pair of real 
vectors v and w obeying + = 1, with the correspondence being given by 

n±i = ^ [^{vx + iwx) - i{vy + iwy)] 
no = Vz + iwz . 

Using Eqs. (4.21) and (5.8) we find that 





1 


= 


An 




3 , 






An 




' ( 


2 = 


207r ^ 



(5.13) 



and hence that 



^i = 7^7^(^P-?2 +^5^ -2P-5gi + g2)|vxw| . (5.14) 



207r V r J 207r V 9 



The nature of the minimum depends on whether the coefficient of |v x w| is positive or 
negative. In the former case, h is minimized when v and w are parallel. By choosing their 
direction to be along the z-axis and then applying a global phase rotation, we can bring the 
solution into the form 



It then follows that 



and that 



riQ — l, n±i = . (5.15) 



3r;2 

r2$*$M = ^sin2^ (5.16) 



with all other (Tjm vanishing. 
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If instead the coefficient is negative, then Ji is minimized when v and w are perpendicular 
and of equal length. Any such solution can be rotated so that Vx = —Wy = — 1/\/2 with all 
other components vanishing. This gives 

ni — 1 , riQ — ri-i — (5.18) 



and 



r2$*$/^ = _(l-cos^)2. (5.19) 



The nonzero cxj^ are 



1 1 

c^00 = ^=, criO = --r^, '^20 = — 7= ■ (5.20) 

V47r 4V7r 4v57r 



Both solutions are axially symmetric; the former is manifestly invariant under a rotation 
about the 2;-axis, while the latter is invariant if the rotation is supplemented by a global 
gauge transformation. 

Case iv: ? = 3, J — 1 

For larger g, the minima of l\ depend on the actual values of the integrals p and gj, 
which we can only determine numerically. However, if A\/g^ is sufficiently large, the first 
term in Eq. (5.7) is dominant and the dependence on the qj can be ignored to leading order. 
In fact, one only has to minimize 



Cjm I 



f (5.21) 
= -J / d<f>d9 sin 9 r"^ K*^)^ • 

The integral in the second line is a sum of integrals of products of four vector harmonics. 
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Using the explicit expressions^ given in Eq. (4.8), we obtain 



^ AMiMiMaMi'^MinM^riM^riM^ (5.22) 

Mi,M2,M3,Mi 



where 



. _^ [(2g-l)!]2(2g + Mi + Af2-2)!(2g-Mi-M2-2)! 

^Mi,M2,M3,M4 — (>{Mi+M2),{M3+M4) / . = ■ 

47r(4g - 3)!y'n,ti(? + Mj - l)!(g - Mj - 1)! 

(5.23) 

Our problem has now been reduced to the minimization of a quartic polynomial in 
2 J + 1 = 5 complex variables. Even after using the rotational and phase freedom to fix some 

of these, one is still left with a rather formidable task. We therefore used Mathematica to 
search for minima, finding a solution that can be rotated into the form 

1 1 
^ ' "-±1 = 0, n±2 = ±- . (5.24) 

Prom this one finds that the nonzero ajm are 



V5 



1 V70 
^ 24A ' ^ 336A 



(5.25) 



This solution has no continuous rotational symmetry, although it is invariant under the group 
of finite rotations that leave the tetrahedron invariant. This is illustrated by Fig. 1, where 
we present a three-dimensional plot of r^^* as a function of angle. Note that vanishes 

at the center of each of the faces of the deformed tetrahedron in this figure, in agreement 
with our previous remark that it should have 2{q — 1) zeros. 



4 Although these expressions for the vector harmonics are gauge-dependent, the result for S is gauge- 
independent. 
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Including the effects of the terms involving the qj shifts the location of the minimum 
of Ji. One can verify that (in contrast with the q = 2 case) the contain terms that are 
linear in the deviations of the um from the values given above. As a result, the full solution 
for the hm changes continuously as X/g^ is varied. 

Since the minimum found here was obtained by numerical methods, we do not have 
an analytic proof that it is in fact the global minimum (although we are fairly confident 
that it is.) However, we can demonstrate unambiguously that the global minimum is not 
axially symmetric. To do this, we note first that any configuration with all but one of the 
Hm equal to zero is invariant under rotations about the z-axis (possibly supplemented by a 
gauge transformation). By evaluating the ajm with m 7^ 0, it is easy to show that these are 
the only configurations with this symmetry. Explicit calculations for the five configurations 
of this form shows that they all give higher values for E than does the configuration of 
(5.24). Hence, the global minimum cannot be achieved by a configuration that is axially 
symmetric about the 2;-axis; the rotational symmetry of the theory then extends this result 
to an arbitrary axis of rotation. 

Larger charges 

We have applied the methods used for the q — 3 case to higher charges also. For q — A 
(i.e., J — 3), the lowest minimum we find for E has 

n±2 = ±^ (5.26) 

with all other um vanishing. The nonzero ajm are 

1 



(700 = 



47r 



70 



13 7V13 

(760 = - r-^r, r- ' ^6,±4 = 



5727? ' 572v^ ■ 

A three-dimensional plot of r^$*$'^ for this solution is shown in Fig. 2. As suggested by the 
plot, this solution is invariant under the discrete rotational symmetries of the cube. has 
a zero on each face of this roughly cubic shape. 
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As wc go to higher values of g, the solutions develop more small-scale structure, while 
at the same time appearing more symmetric when viewed on a large scale. (Note that a 
discrete polyhedral symmetry such as that exhibited by the g = 3 and q = A solutions is 
impossible for most values of q.) What we see happening is that there is a tendency for the 
2{q — 1) zeros of to be distributed as evenly as possible over the unit two-sphere. Lying 
between these zeros are maxima of This behavior can be seen, for example, in Fig. 3, 

where we show a solution with q — 12. 



The deviation hfj,^, of the metric from the Reissner-Nordstrom solution is determined to 
leading order by the linearized Einstein equation 



Here denotes the terms in the Einstein tensor that are linear in h^i, while t^,^ is the 

leading correction to the energy-momentum tensor. 

In doing this calculation, we continue to restrict ourselves to the case where Gw? /e^ is 
very small; it was this assumption that allowed us to decouple the determination of SF^y from 
that of h^y. For our perturbative scheme to be valid, the horizon radius of the unperturbed 
solution must be close to rgr, and hence must be of order . When this is the case, the 
difference between the Reissner-Nordstrom and Schwartzschild metrics with the same value 
for M is never greater than order Gm^ / anywhere outside the horizon. Hence, in solving 
for h^jti, we can approximate the metric by the corresponding Schwarzschild metric. We will 
consider only the case g > 1, > 0, so that we can use the results for the VF-field and the 
electromagnetic perturbations that were obtained in Sees. 4 and 5. 

The first step is to calculate t^u. The full energy-momentum tensor can be written as 



6. Perturbation of the Metric 



5Gau = -STrGt 



(6.1) 




nEM I rpW 
/if /ttf 



(6.2) 
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where 

is the purely electromagnetic part and 



(6.3) 



(6.4) 



The 0{a^) corrections to T^/^ are the sum of a part linear in ?)F^i, and a part linear in h^J^l, 
that arises from the corrections to the metric in Eq. (6.3); because the latter is suppressed by 
an additional factor of Gm^a^, we can ignore it here. The dominant part of TJ^, which is also 
O(a^), is obtained by substituting the unperturbed metric and field strength into Eq. (6.4). 
Using Eq. (1.5) for the only nonzero component of the unperturbed electromagetic field 
strength, together with the results of Sec. 4, we find that the nonzero components of t^i, can 
be written as 

tu = B {Kf, + [Biff + {m' - ^) f] K,y} 

trr = ~ {Kfi + [-Biff + (m' - ^) f] Ky,} 



tab = 9ab l-f^Fl - TToJ J^W 



(6.5) 



2r2- 



Kf2 



where 



Kw 



Kfi 



mr 



jm 



YjmiO,4>) 



er 



^1 ^jm 



jm 



^F2- -^Yl F^^daYjm - -^-^-^^''jfJjmdaYjm 



jm 
qga^ 



(6.6) 



jm jm 

(In the last two lines we have used Eqs. (4.30) and (4.32) to relate the perturbations of the 
field strengths to those of W^; it should be recalled that = 0.) 
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The next step is to expand the components of hf^^, in terms of spherical harmonics. The 
space-space components of /i^j^ can be decomposed into a spin-0 field and a spin-2 field, 
the time-time component corresponds to a spin-0 field, and time-space components can be 
chosen to vanish because the solution is static. Thus, there are potentially seven functions of 
r entering this expansion for each value of j and m. However, examination of the parity of 
and t^j^ shows that it is sufficient to consider only those terms corresponding to perturbations 
of parity (— l)-^- (In the terminology of Ref. 12, these are polar perturbations.) In general, 
this leaves only five modes for each value of j and m: two, with I — j, for the spin-0 fields, 
and three, with I = j — 2, j and j ' + 2, for the spin-2 field. Hence, there are five radial 
functions, which we define by 



jm 



jm 



jm 

hra = Y Hl'^{r)VaYjmK'y, ' 



jm 

where denotes the generally covariant derivative with respect to the unperturbed metric. 
For j — 1 there is no mode with I — j — 2, and so there should be only four radial functions. 
Indeed, the identity 

VaVbYiM{0, (t>) = -^YiM{0, <P) (6.8) 

shows that the Hl'^ is redundant and can be set equal to zero. Similarly, there should be 

only three radial functions for j — 0, and the fact that Yqq is a constant allows us to set 
i/oo ^ ^00 ^ 0. 

Further simplification can be achieved by utilizing the freedom to perform coordinate 
transformations, which change the metric by an amount 

SGQfiu = V^ei. + V^e^ . (6.9) 
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Writing 



et^O 

jm 
jm 



(6.10) 



we find that 



dcH-r = -B'fr 
5GHi'^=2Bf'r+B'fr 

SGHr = ^-^fr (6.11) 
scHr = 2/r 

soHt = /r + fT - Ifr ■ 

For j > 2, we choose Z^"" so that Hi"" = and then so that hI"^ = 0. For j = 1, we 
choose fl"^ so that HI"^ = and then choose /g"' so that ifg"'' = 0. FinaUy, for j = 
we choose /{'^ (the only coordinate freedom available) to set H^^ = 0. 

We next note that tab is proportional to Qab, even though it could in principle have also 
contained terms proportional to VaV&FjM(^, <P)- The absence of such terms implies that 

sin2 9 5Gee - SG^^ = ^ J] (//^ - Hf) (sin^ 9 - V^V^) Yjm = (6.12) 

jm 

from which it follows that 

HJrn^jjpn^ j>2. (6.13) 



This leaves only two independent radial functions for each value of j and m and gives 



29 



us a metric of the form 



gtt = B{r) 



9rr 



Bir) 



j=0 m 



1 + i/2°° W>oo + E E Htir)Yjm{0, 4>) 

j=l m 



(6.14) 



5*06 = 9ab 



i+EE^rw--(^,0) 

j=l m 



All of the must vanish as r — > 00. At the horizon, the tt and ab components of the 

metric, as well as its determinant, are nonsingular. Wc require the same of the perturbed 
metric, and hence require that BH^^, iJg^ — iJj'^, and iJg"^ all be nonsingular at r = th- We 
will see that the field equations place further restrictions on the behavior near the horizon. 

Differential equations for the various H^J^ are obtained by expanding both sides of the 
linearized Einstein equation (6.1) in terms of spherical harmonics. Because of the Bianchi 
identity, as well as the symmetry of the problem, many components of the resulting equations 
are redundant. In particular, if we write 



5Gu = Y.5GlT{r)Y^m{e,^) 



5Gra 



jm 



jm 



(6.15) 



jm 



then it is sufficient to calculate 
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B 



B 



[(1 - tB' - B){H^^ + Hf) - - TB{Hf)'] 
{rBHf)' + 0{Grr?le^) 



r^B 
1 



r^B' 



rBiHi'^y -lrB+-^] {HPY + ^ ^^ H^ - Hi"^) 



{Hf^y - {Hi"")' 



y I (j-l)(j + 2) 
' 2 

J >1- 



(6.16) 

In the second equality for SG^^ we have used the fact that, with the approximations we are 
making, 

1 - rB'{r) - B(r) = 0{Gm^/e^) (6.17) 

everywhere outside the horizon; for a Schwartzschild metric the left hand side of this equation 
would vanish identically. 

We start with the j = modes, for which we can use Eq. (4.22). To leading order, the 
tt component of Eq. (6.1) leads to 



(6.18) 



A second equation is obtained by multiplying the rr equation by B{r) and the tt equation 
by 1/B{r), and then adding these to give 



47r TTir 



(6.19) 



These two equations can be immediately integrated, using the boundary condition that 
Hf^{oo) = i/f (oo) = 0, to give 



S(r)//f(r) = 



2x/A^G6M 2x/4^Ga2 



mr 



ds 



B 



{s){f'{s)f + {m'-^)f{s)] (6.20) 
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B{r)H^°{r) = — / ds 



mr 



B{s) 



2B(r)r 



(6.21) 

where SM is an arbitrary constant that may be interpreted as a shift of the black hole mass. 
Because /(r) falls exponentially fast at large distance, so do the integrals on the right hand 
sides of these two equations. 

Turning now to the j ^ modes, we find that the ra components of Eq. (6.1) give 



BJ^'; 



J>1 



(6.22) 



while the rr component leads to 



= 8.GaV,. {1 [-Biff + {m^ - ^) f] + ^^,} 



(6.23) 



It was noted above that the nonsingularity of gu required that BH^ be regular at r = rj^; 
this leaves the possibihty that might be singular there. This can be ruled out by 

multiplying Eq. (6.22) by r and then adding the result to Eq. (6.23). The resulting equation 
can be solved to express in terms of H^^ and [H^^)' , thus showing that it is regular 
at the horizon. Hence, (iJ^"^)' is finite at ru- Using this fact in Eq. (6.22), we find that 



Hi"\rH)^Q, J>1 



(6.24) 



For j = 1, the fact that i/g only enters Eqs. (6.22) and (6.23) through its derivative 
simplifies matters considerably. By using Eq. (6.23) to solve for {H^^)' and then substituting 
into Eq. (6.22), we obtain a first-order equation involving only Hf^. This can be easily 
integrated, and the result then used to obtain H^^. The two constants of integration are 
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fixed by the boundary conditions that Hl^^rn) = H^^{oo) = 0. The result is that 



(6.25) 



+- 



m 



H^"^^ J ds l'^Hl'^{s) + 



f (5i?(.)^((.)H-^i(.)) 



+- 



ms 



sB'(s) 

5(.)(/W+(-^-f)/ 



(6.26) 



To obtain the asymptotic behavior of these expressions, we recall that /(r) vanishes expo- 
nentially fast at large r, while J-'i ~ 1/r. This behavior is just sufficient to guarantee the 
convergence of the integral in Eq. (6.25), with the result that ~ Inserting this 

into Eq. (6.26) and again using the asymptotic behavior of Tj, we find that H^"^ has the 
same asymptotic behavior, and in fact that the difference HI'^ — ~ 

In solving for the modes with j > 2, it is useful to define 



(6.27) 



and to then rewrite Eq. (6.22) as 



rrjm _ _jB 

^1 " B' 



3'm 



19'^ jm \ 



j{j + l)J r 



(6.28) 



Substitution of this equation and its derivative into Eq. (6.23) leads to a second order equa- 
tion involving only Tjm'- 



r^BT'f + 2{r'B)'T'^ - {j - l)(j + 2)7}-^ = 



jm 



m 



[-B{f'f + m^f] . (6.29) 



(In obtaining this result, both Eqs. (4.33) and (6.17) have been used, and terms of higher 
order in Cm? /e^ dropped.) 
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In the absence of the source term there are no nontrivial solutions for Tjm that are regular 
at both the horizon and spatial infinity. (This can be readily seen by examining the explicit 
solutions given in the appendix.) Hence, the inhomeogeneous Eq. (6.29) can be solved by 
Green's function methods similar to those used for Eq. (4.33). Let h'j{r) and /i^(r) be the 
solutions of the homogeneous equation that are regular at r ~ rjj and r — oo, respectively. 
At large r these behave asymptotically as r^~^ and r~^^'^'^\ If they are normahzed so that 
r~^^~^^h~{r) and r^'^'^hj'{r) both tend to unity as r — > oo, then the desired Green's function 
is 

Cf Vr,r') = -TTT^ \e(r - r')ht (r)hT (r') + e(r' - r)hj (r)h+ (r')] (6.30) 

J ^7 ~\~ 1 J J J J 

and the only regular solution for Tjm is 

Tjmir) = J dr'Gf\ry){rfB{r') {-B{r')[f{r')]' + m'[f{r')]'} . (6.31) 

can be obtained immediately from this equation together with Eqs. (4.34) and (6.28). 

At large r the source term in Eq. (6.29) is exponentially small, so the first term in the 
Green's function dominates Eq. (6.31); thus 

C- 

Tjm ~ (^jm r ^ OO (6.32) 



where 

Gj — 



(2j + 



'^Jdrr'BhjiBiff-m'f] . (6.33) 

By substituting this result back into Eq. (6.28) and using Eq. (4.36) we obtain the asymptotic 
behavior 

Tjjm SnGa^ajm Aj 

Hi ~ T^T, r ^ oo 6.34 

1 {2j + l)mrHr3+^ 



where 



Aj -Jdr\ (2(j + 2yBh- [B{ff - m'f] + qgj^-^ I . (6.35) 
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The resulting large distance behavior of gu corresponds to that one would obtain from a mass 
distribution with a 2-' -pole moment whose components are equal to the ajm times quantities 
of order p'^r-'jf^ / ~ /9^/e^m-^~^. 



7. Concluding Remarks 



In this paper we have exhibited black hole solutions with fields on the horizon that, 
contrary to common expectation, are not spherically symmetric. As the magnetic charge 
increases, the structure of these black holes becomes more detailed, with higher multipole 
components appearing in the long-range electromagnetic and gravitational fields. At the 
same time, the lower multipole moments decrease in magnitude, with the result that the 
solutions do begin to approach the expected spherical symmetry, if only in an averaged 
sense. 

Although these solutions display some unusual, and perhaps unexpected, properties, 
there are rigorous general results on black holes that they must obey. We consider here two 
of these, the zeroth and second laws of black hole dynamics, that are concerned with the 
properties of the black hole horizon. At constant t, the horizon can be described as the 
two-dimensional surface 



where rn is the horizon radius of the unperturbed metric. To leading order, the vanishing 
of gtt on this surface gives 



jm 

Eq. (6.24) implies that the terms in the sum with j > 1 all vanish, so that A is a constant 
independent of angle. Using the solution in Eq. (6.21) for H^'^ (with dM = 0), we find 



r = rH + A{9, (p) 





(7.2) 



CO 



2Ga'^ 




ds [Biff + (m^ - I) f 



(7.3) 
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where the second equahty follows from the eigenvalue and normalization Eqs. (4.10) and 
(4.11), as well as Eq. (6.17). Because the integral in the first line is itself of order P'^, our 
result for A is proportional to rather than the one might have expected. Since the 
the higher order corrections to FF^ and SFf^i, could also shift the horizon by a distance of 
order our result is really that A vanishes to leading order. 

The zeroth law of black hole dynamics states that the surface gravity — which corre- 
sponds quantum mechanically to the black hole temperature — is constant over the horizon. 
For a stationary black hole the surface gravity k. is given by [13] 

«^' = -^(V^X^)(V^X.) (7.4) 

where the right-hand side is to be evaluated on the horizon and Xn is a Killing vector that 
is orthogonal to the horizon. If the metric is actually static, and is written in a manifestly 
^-independent form with vanishing time-space components gtj — 0, then the only nonzero 
component of this KiUing vector is xt — Qtt and Eq. (7.4) reduces to 

^''^-\9''9''{digumgu). (7.5) 

Expanding this equation to first order in h^jj and taking into account the fact that the 
horizon has been shifted by an amount A gives 

B' i?"A B' 
« = Y + drhtt + — {B~^htt - Bhrr) 

B' B"A l^.„„,w„, B'. .„nn „nnx ^^'^^ 



2 2 2 

jm 

where all quantities are to be evaluated at r = rn- Because of Eq. (6.24), all terms in the 
sum with j > I vanish. Hence, k is independent of angle and therefore constant over the 
horizon, as required. 

The second law of black hole dynamics is the statement that the area of a black hole 
horizon never decreases. Let us apply this to the case of an unstable Reissner-Nordstrom 
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black hole that is perturbed and eventually evolves into a static black hole with hair. By 
making the perturbation sufficiently weak, we can arrange that the mass decrease from 
radiation be negligible, so that the final state will be a black hole with hair that has the 
same value for M as the original Reissner- Nordstrom solution. (Indeed, the area law can be 
used to place an upper limit on the mass loss from radiation.) The area of its horizon is 



where the integration is over the surface r = r^y + A. The angular components hab of the 
metric perturbation only have terms involving involving spherical harmonics with j > 1. The 
contributions linear in these vanish after the integration over angles, and so for calculating 
the leading correction to the area we can replace gab by the unperturbed metric gab and 
obtain 



Prom this, together with Eq. (7.3), we see that the second law is verified, to leading order. 

The methods we have used to construct our solutions have limited us to the case where 
the horizon radius is close to the critical radius for the instability of the Reissner- Nordstrom 
solution. However, there seems to be every reason to expect that the solutions with smaller 
horizons will display similar behavior. The construction and study of such solutions, which 
include new extremal black holes, remains an open problem. 




(7.7) 



A = 47r(r|f + 2riy A) . 



(7.8) 
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APPENDIX 



The Green's functions used to solve Eqs. (4.33) and (6.29) were constructed from the 
solutions of the corresponding homogeneous solutions. In the former case, the solutions obey 



B{T)g"{T) + B'{T)g'{T) - ^^^^(r) = . 



(A.1) 



As we explain in Sec. 6, in the approximation to which we are working B{r) may be replaced 
by the Schwarzschild metric function l — ru/r. It is convenient to define a variable x — r/ru 
and rewrite Eq. (A.l) as 



This has a polynomial solution of the form 



(A.2) 



(A.3) 



where P^^^ {z) is a Jacobi polynomial [14], and the normalization has been chosen so that 



^ ' 9j (^) tends to unity as a; oo. Its value at x = 1 (i.e., the horizon) is 



^7(1) = 



(J-1)!(J + 1)! 
(2j)! 



(A.4) 



To find the other independent solution, we first note that any two solutions g and / of 
Eq. (A.l) obey 

B{r)[g'{r)f{r)-g{r)nr)\^c (A.5) 
where c is a constant. In particular, if we write 



g+[x) = 5f • (x) In ( 1 - - ) + k{x 



X 



(A.6) 



then k{x) obeys 



X — 1 



X 



_dk ^dg- 
dx dx 



+ 



{9, 



c . 



x^ 



(A.7) 



From the fact that g^ (x) is equal to x^ times a polynomial of order (j — 1), it follows that 
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this last equation has a solution for k{x) as a polynomial of order j. 
Explicit forms for low values of j are 



9i {x) = x'^ 



^2"(^)=^'(^-i) (A.8) 



,2 / .2 4 2 



and 



g+{x)^-3x^\n(^l-^^ -Sx-^ 

g+{x) = -80^2 (^x - ^ In (^1-^^ - SOx^ + 20x + Y (A.9) 

, , , o / o 4 2\ / 1\ o 2625 o 35 

g+{x) = -1575x2 ( 3^ ^ 5 ) ( ^ ~ X ) ~ ^^''^^ + ~ -^^^^ ~ T 

where the g^{x) have been normalized so that x^gj~{x) tends to unity as x — > 00. 

With the same approximation for B, the homogeneous equation corresponding to Eq. (6.29) 
becomes 

This has a polynomial solution of the form 



(x^ - x)-^ + (4x - 2)'=- - (j - l)(j + 2)ft = 0. (A.IO) 



h-(x) = ^^^^^l^Pj]:?(2x-l) (A.11) 
whose value at the horizon is 

"7(1) = '-^1^1^. (A.12) 

By methods similar to those used to find gj'{x), one finds that the solution that is regular 
as X —> 00 is of the form 

A;W = /.-M1.(i-1)+^ (A.13) 

where £{x) is a jth order polynomial. 
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Explicit forms for low values of j are 




(A.14) 



/ig (X) — — X + 



and 



ht(x) = 60(2x - 1) In (1 - -)+ 120 - , , 

\ X J x{x - 1) 




(A.15) 



h+{x) = 35 1 12(5x2 - 5x + 1) In ( 1 - i ) + 60x - 30 + 4^ 




where the normalization conventions are analogous to those used for g- (x) . 
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FIGURE CAPTIONS 

1) A three-dimensional spherical plot of the quantity r^$*$^ for the g = 3 solution 
described in Eq. (5.24). If we denote the spherical coordinates of a point as {R, 6, 4>), 
then this plot shows the surface R{9,(p) = r2$*(^,0)$/^(^,0). Note that while 9 
and (f) represent the corresponding spatial coordinates, R is unrelated to any physical 
spacetime coordinate. Note the tetrahedral symmetry of the surface. 

2) A spherical plot, similar to that shown in Fig. 1, for the q — A solution of Eq. (5.26). 
The cubical symmetry of the solution is apparent. 

3) (a) A spherical plot, similar to those in Figs. 1 and 2, for a solution with q = 12. 
There is no apparent symmetry. 

(b) Another presentation of the same solution. The value of the function r^$* on 
the unit hemisphere (0 < < tt) is represented by one of 16 gray levels, with black 
being the minimum (zero) and white being the maximum. One can see the fairly even 
distribution of the zeros. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/gr-qc/9503035vl 
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